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Abstract. The scale invariant properties of wave functions in finite samples of one 
dimensional random systems with correlated disorder are analyzed. The random 
dimer model and its generalizations are considered and the wave functions are 
compared. Generalized entropic localization lengths are introduced in order to 
characterize the states and compared with their behavior for exponential localization. 
An acceptable agreement is obtained, however, the exponential form seems to be an 
oversimplification in the presence of correlated disorder. According to our analysis 
in the case of the random dimer model and the two new models the possibility of 
power-law localization cannot be ruled out. 



PACS numbers: 71.23.An, 72.15.Rn, 73.20.Jc 
1. Introduction 

In a previous publication [JjJ we have introduced a new form of information length 
in order to characterize the shape of wave functions in finite one dimensional (Id) 
disordered systems. Using that definition we succeeded to show that the states in the 
Id Anderson model with uncorrelated, onsite disorder do have, apart from oscillations, 
an overall exponential shape. Such an exponential decay has been found for practically 
any strength of disorder even in the case when the localization length exceeded the 
size of the system substantially. 

The scaling properties of one-particle states in the presence of uncorrelated 
disorder in Id and quasi-ld has been studied extensively both numerically and 
analytically |[ ||, |). The similar problem of the more realistic case of correlated 
disorder has been recently considered in [0, fj\. In this paper we wish to present a 
scale invariant study on a wider family of correlated disorder in Id and at the same 
time show how generalized localization lengths may help to analyze the properties of 
the one-particle eigenstates. 

To be more specific the eigenvalue problem of an electron in a Id disordered 
potential can be given as 

E Cm — £"m C-m ~\~ ^rr^m+l ('m-\-l ~\~ Vm — l,m (1) 

where c m is the amplitude for the electron to be on site m and E is the energy 
eigenvalue. In the most simple case studied in |l| the onsite potentials e m are chosen 
randomly from a box distribution of width W centered around the origin, and the 
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off-diagonal hopping integrals are kept constant, V mm +i = V m -\,m = Vo- The latter 
condition enables us to fix the unit of the energy scale Vq = 1. For this model there 
are rigorous results || affirming complete exponential localization for any strength of 
disorder in infinite systems and in Q we have proved numerically the above statement 
holds in finite systems, as well. 

The effect of correlated disorder introduced in (|l]) has attracted much attention 
recently. These correlations maybe originated from interactions of electrons with 
lattice vibrations or e.g. as a more realistic representation of disorder incorporating 
the presence of chemical bonding. The first models that included such correlations 
were therefore based on random binary models §, |§ £| H Q. These studies 
revealed the possibility that disorder correlations may increase and as well as decrease 
the localization length substantially. 

In a special family of random binary alloys correlation was introduced by assigning 
the same energy level ea or eb to pairs of sites. This is called the random binary 
dimer model (RBDM). This model was first introduced in and it has been shown 
that under certain conditions there are special E c values where the state is delocalized, 
transparent and the number of extended states around E c is proportional to the square 
root of the length of a finite sample. The localization length diverges at these energies 
H and this is reflected in the conduction properties of finite samples jL6). In a 
recent paper JtJ Izrailev et al. have studied the scaling properties of the eigenstates 
in the RBDM and succeeded to show that the states approaching E c are described 
similarly as in the case of uncorrelated disorder jl], ||. 

There are other similar models formulated in the same spirit as the RBDM. A 
continuous Kronig-Penney type random dimer model |l7t exhibits an infinite number 
of resonances with zero reflection constant. The existence of similar special states in 
a quasiperiodic dimer model has also been found in [ jl8| . 

Exponential localization, although with an enhanced localization length has been 
seen in another model where the e m energies are drawn from a box distribution but 
they are repeated for L consecutive sites [ p!9[ . Here it is possible to vary L, however, 
no special energy with complete derealization is present. 

In this paper we present numerical results of two generalized versions (A, B) of 
the RBDM that are related to both the model of finite correlation length in [1^] and 
the original model of Dunlap, Kundu and Phillips [Q. Our results are compared to 
the ones obtained by Izrailev et al. 0. 

In Model A the onsite energies are drawn from a box distribution and assigned to 
two consecutive sites at the same time: it may be called general random dimer model 
(GRDM). This model is intermediate between the original Anderson model and the 
RBDM, it is in fact the special case of L — 2 of the model studied in As it has been 
shown fH| within such models the energy band does not contain any special energies 
where complete derealization may occur, however, correlations change localization in 
a direction similar to the RBDM. 

Model B |2(| on the other hand contains disorder in both the diagonal and the 
off-diagonal part of 

£m = — Vo(a m ,m+l + Q; m ,m-l) 
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Ki,m±i = Vby 1 + ot 2 mm±1 - 2a m . m ±iCOs5 (2) 
where the quantities a m ,m±i are chosen from a box distribution centered around the 
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origin with width W < 2. This model is obtained considering the coupling of electrons 
to the vibrations of the underlying lattice represented by the random variables a m) m±i 
that introduce a correlated disorder in both the onsite and the off-diagonal matrix 
elements (see Q for the details). The correlation is perfect if in ^) G = 7. The 
RBDM can be considered as a simplified version of this model. The special energies 
where derealization occurs are at E c — 2Vq cos S. We chose Vq as the unit of energy 
here and varied the energy in the vicinity of E c for different values of 8. 

The solution of the Schrodinger equation ([!]) using appropriate initial conditions 
c a = c\ = 1 is obtained by numerically iterating for a system of N = 10 4 sites and the 
localization properties of the eigenstates are calculated using the charge distribution 
Qm = l c m| 2 - Averaging is performed over M = 1000 samples. 



2. Shape analysis 

The shape of the charge distribution Q m may be characterized using the inverse 
participation number (IPN), D, and the Shannon-entropy, H |Q] 

D- l =Y J Q 2 m and H = -J2 QmlnQm. (3) 

rn m 

Both parameters D and exp(H) give the number of sites effectively populated by the 
state. Therefore a state extending over the whole system would have D = N and 
H = InN. This means that one may introduce two parameters [Q 

Pi = cxp(i7 - H n f) and 2 = D/D ref , (4) 

where normalization with respect to both the system size and the case of the absence of 
disorder has been performed. The latter is achieved by evaluating the values D re f and 
H re f for the Bloch-wave solution of with e m = and V m , m ±i = Vq Jj- Overbar 
means average over the samples at fixed energy and/or disorder. It is clear that both 
of these quantities change from 1 to as disorder increases from to 00 therefore in 
they have been used as scaling functions. 

Another way of expressing pi and Pi (^) can be obtained using our previous 
definitions [Q. There we have used a size independent form and pointed out its 
relevance in a scale independent shape analysis of the states. In we have calculated 
the spatial filling factor, q, and the structural entropy, S of individual eigenstates as 

q = D/N and S str = H - In D. (5) 

These quantities obey the inequalities: < q < 1 and < S s t r < — In q. In the absence 
of disorder the solution of (|l]) is a plane wave for which q° — 2/3 and S*° tr = In 3 — 1. 
These quantities can be related to the reference values in (|4|) as q° = D re f/N and 
Sg tr = H re f —\nD re f. Note that q° and S° tr are independent of the system size while 
D re f and H re f are not. Using q and S s t r we may rewrite equations (^) in the form 

Pi = q exp(S str ) / Po and p 2 = q/q°, (6) 

where po = q° exp(S® tr ) ~ 0.7357. We have to note that the quantities q°, Sg tr and 
Po that appear in (^) are obtained naturally from the solution of ([!]) for uncorrelated 
disorder which is a plane wave modulated by an exponentially decaying envelope 
represented as a product: c m ~ exp(— |m — mo \/Qsm(km + 5) g. The very same 
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conclusion was drawn using a completely different method by Fyodorov and Mirlin j5| 
based on results of quasi-ld systems and that of strictly Id j2jj . 

The main advantage of this reformulation is the application of the shape analysis 
proposed originally in p2[ . We have shown in Q and in other publications [^3[ p4fl 
that our method is applicable for eigenstates composed as a product of several simple 
forms, e.g. an oscillating planewave and an envelope characterized by the scale £. In 
p3[ our method unambiguously showed the existence of power-law delocalized states 
near the mobiliy edge of a Id quasiperiodic system, where both the bulk and the tail 
of the envelope played an equal role in the analysis of the wave functions. 

It has been shown in [p2[ that parameters q and S str can be calculated for ideal 
charge distributions analytically and that the relation S str (q) is directly connected to 
the shape of the distribution. Therefore the properties of a large set of numerically 
obtained wave functions are compared to ideal curves in the parameter space (q, S s t r ) 
especially when some control parameters of the system are varied, e.g. system size, 
strength of disorder or energy. Similar relation may follow between f3\ and /fe- For 
example in the case of exponential localization we obtain 

ft( Z ) = *LieJl * ) ^ 

z exp(z) \ exp(z) — 1 J 



and 



z \ exp z + 1 

where z = N/£ with localization length £. [3\{z) and 0%{z) are monotonous functions 
of z, hence the relation 0i — > /?2 exists and is also directly connected to the shape of 
the charge distribution. Note that similar analytic (3\{z) and ^{z) functions can be 
calculated for any other type of form functions, e.g. for power-law decay. 

3. Results and discussion 

In Izrailev et al. fitted a very simple analytic form for the relation /3i — » /?2 in the 
case of the RBDM 

with c w 0.5488. We will show that this relation is a good approximation, however, 
it fails to describe the states in both the localized and delocalized limits. We have to 
note that a similarly simple scaling relation || for the case of uncorrelated disorder 
has been exhaustively studied in [^| and shown the limitations of it. 
In terms of q and S s t r relation (^) reads as 

S str = -Hc+(l-c)q] (10) 

where S s t r = S s t r — S® tr and q — /?2 = q/q°- First of all in the most extended limit, 
q — * 1, S s tr(q) ~ (1 — q)/2 should hold [^J. In contrast, according to ( |l0| ) we get 
S str (g)«(l-c)(l-g). 

Secondly we plotted the results of our calculations together with the data obtained 
from [Q in Figure [ll The analytical forms by equations (solid curve) and also the 
empirical relation (^) (dotted curve) are shown, as well. It is clear that (^) is indeed 
a good approximation, however, the tendency is somewhat closer to equations (0,||) 
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Figure 1. Interrelation Pz(ft\) of the generalized entropic lengths. Solid symbols arc 
from J?], the crosses (DKPM) and open squares (GRDM) are results of the present 
calculation. The continuous curves stand for analytical relations. In the inset the 
difference A/3 = P2 — Pi is given as a function of /3i. See details in the text. 



that shows exponential localization on all length scales. The third relation (dashed 
curve) is the one derived assuming an envelope of the form c m ~ m~ 3 instead of an 
exponential form. The inset shows the deviation A/3 = /3 2 — /3i as a function of j3\. 
In the inset we see again that relation (^) is an acceptable fit to the data from j7j, 
however, in this figure it is very hard to check its accuracy especially for the localized 
and delocalized limits. In Figure [j] we have plotted our results for Models A and 
B, as well. For Model A we have varied the width of the box distribution between 
W = KTVo up to W = 10 4 Vb. In the case of Model B the parameters W = 1.9V 
(here W is limited to < W < 2) and S = 0°, 30°, 60°, 90° have been used. We can 
see that, at least in terms of the (3i — > 02 relation, both of these models are not very 
much different from the behavior of RBDM. 

In order to investigate the similarities and differences between the RBDM and the 
models studied here (Model A, the GRDM and Model B, the DKPM) it is even more 
apparent to plot S str as a function of In q. In Figure ^ the localized region q — > 
is clearly not described by the relation ([l0]) presented with a dotted line. It is also 
true that the states in neither models show clear exponential localization in the ideal 
form depicted by the continuous solid curve as it is already suggested in the previous 
paragraph. However, the states in Models A and B studied here are closer to it. 
This means that in the strong localization limit, q — > 0, the states in the RBDM have 
definitely more complex structure than that of Models A and B. It is interesting to note 
that similar deviation has been obtained for the case of weak uncorrelated disorder 
in pH |. However, in contrast to pi] ] our results are nonperturbative. Furthermore we 
have to stress that in the present work disorder correlations play an important role 
yielding the above mentioned deviations from conventional exponential localization. 
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Figure 2. The localization diagram for states in the RBDM, GRDM, and the 

DKPM. The relation (nut) ( ) clearly deviates from the numerical results for 

q — > 0. The symbols are like in Figure hi In the inset the deviation of S a tr is given 
from its expected universal form (1 — q)/2 when g — » 1. The dotted line is wrong in 
this limit, as well. 



In Figure || we have also plotted the S s t r {q) relation for power-law localization 
with different exponents. We observe that the RBDM is well described by an overall 
shape: c m ~ m~ 3 (dashed line) while the GRDM and the DKPM are better described 
with c m ~ m~ 6 (dashed-dotted line). Apparently this is in contradiction with 
analytical expressions of the Lyapunov exponent (inverse localization length) which 
goes as 7(i?) ~ (E — E c ) 2 around the special energies E c |q, O, however, 7 should 



vanish for the case of power-law localization [|23| 24 . A possible resolution to this 
problem is already outlined in Section e.g. an exponential decay superimposed by 
some kind of rapidly varying substructure can easily provide a shift from the curve 
corresponding to the exponential decay to power-law decay as it is seen in Figure ||. 
According to Figure ^| the relation (10) is a wrong approximation especially for the 
strong localization limit. Moreover, in the inset of Figure^ we see deviations for the 
delocalized limit, q — > 1, as well. 



4. Conclusions 

We have performed a shape analysis of wave functions obtained in several one 
dimensional random models with correlated disorder. We have introduced a new 
definition for the generalized localization lengths based on the inverse participation 
number and the Shannon-entropy. We have applied the shape analysis introduced in 

It has been shown that the localization properties of the states in the RBDM are 
described by (||) only approximately. On the other hand equation ( |l0| ) shows a wrong 
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behavior for q — > 1. Large deviations are obtained in the localized limit /3 2 — ► 0. 
We have compared the data from j?J with our simulations on the GRDM (Model A) 
and the DKPM (Model B). The data show a clear deviation from simple exponential 
localization: the average localization properties of the states for the RBDM resemble 
that of a power-law shape c m ~ m -3 and in the case of Models A and B that of a 
power-law shape c m ~ to -6 . 
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